AD-A168 899 7 SEMI-EMPIRICAL LON-LATITUDE IOOSPHERIC MODELCU) AIR 
FORCE GEOPHYSICS LAB HANSCOM 
D N ANDERSON ET ale 1@ OCT 835 AFGL-TR-85-0254 
UNCLASSIFIED F/ 


ee oa 


ave eo 8 


tact 


BO Ry PRE yet se ton 


ea 


aoe 


Seng 


= 
w 
o 


ry 
° 


TTP cee 
i 


rr 
r 
te 


i 


| 


less 


SEE 
N 
© 
%S 
Cs 


is 


BN 
WO 


lle 


ks 


li 


AFGL-TR-85-0254 
ENVIRONMENTAL RESEARCH PAPERS, NO. 933 


AD-A168 899 


A Semi-Empirical, Low-Latitude lonospheric Model 


DAVID N. ANDERSON D i ie. * 
MICHAEL MENDILLO : 
BRUCE HERNITER ELECTE 3 

JUN 2 3 066 


D 


10 October 1985 


Approved for public release; distribution unlimited. 


IONOSPHERIC PHYSICS DIVISION PROJECT 2310 


AIR FORCE GEOPHYSICS LABORATORY 
HANSCOM AFB, MA 01731 


g6 6 23 039 


bees Meee fag te aloe Fight ee) Fae ee ahs EE Teh eR OS Se, ety | at 
Pee, ee ee ea 


eS ergs Sn Fie, 


oar ihe See Sera 


é . 
= 
4, 
¢ 
of 
- 
7, 
‘- 
re 
“ 
‘% 
~ ‘ 
ut 
“This technical report has been reviewed and is approved for publication." 
, gr oe 
a 
t 
FOR THE COMMANDER A ce 
HERBERT C, CARLSON Ss 
Branch Chief ge 
en 
This document has been reviewed by the ESD Public Affairs Office (PA) and is os 
releasable to the National Technical Information Service (NTIS). * 
Qualified requestors may obtain additional copies from the Defense Technical a 
Information Center, Al) others should apply to the National Technical Information 
Service. _ 


ig If your address has changed, or if you wish to be removed from the mailing list, on 
or if the addressee is no longer employed by your organization, please notify ae 
AFGL/DAA, Hanscom AFB, MA 01731. This will assist us in maintaining a current ‘. 
mailing list. . 


f i 
. 


a 
’ 
« 

' 

r 
. 
ry 

Ca 
1, 
Ls 

o 


£E en. 


Cee 


Unclassified 


SECURITY CLASSIFICATION OF THIS PAGE ll ‘J G4 IP 


REPORT DOCUMENTATION PAGE 


1a REPORT SECURITY CLASS!FICATION 1b. RESTRICTIVE MARKINGS 
Unclassified 


ft, 


i ee 


ea | 


2a SECURITY CLASSIFICATION MUTHORITY 3 DISTAIBUTION-AVAILABILITY OF REPORT a 
Approved for public release; nt 
y 
4 PERFORMING ORGANIZATION REPORT NUMBER(S) RING ORGANIZATION REPORT NUMBERS) 
AFGL-TR-85 -0254 = 
ERP, No, 933 = 
6a NAME OF PERFORMING ORGANIZATION Bo OFFICE SYMBOL Ja NAME OF MONITORING ORGANIZATION 
c 


‘ \ ey + Hf applicable? 
Air Force Geophysics ; rs 
Laboratory AS 
6c ADDRESS (City. State and ZIP Codes 


Hanscom AFB 


To AOORESS (City, State and ZIP Code! 


° Massachusetts 017: i 
Ba NAME OF FUNDING. SPONSORING Bo OFFICE SYMBOL 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER es 
ORGANIZATION thf applicable, ae 
Bo ADDRESS City. State and 71 Code! 10 SOURCE OF FUNDING NOS 
PROGRAM PROJECT WORK UNIT 
ELEMENT NO NO NO 
TT TITLE Include Security Classificonony A Semi -Empirical 61102} 2310 ; 01 2g 
, Low -Latitude Ionospheric Model “S 
, 12 PERSONAL AUTHORIS) c “ A 
David No Anderson; Michael Mendillo; © Bruce Herniter pare 
13a. TYPE OF REPORT 13b TIME COVERED 14 DATE OF REPORT Yr Mo. Days - 
Scientific, Interim, FROM __. fo LTO ge ek 19835 October 10 “y 
16 SUPPLEMENTARY NOTATION 
Department of Astronomy, Boston Universitv, Boston, MA 02215 oy 
. ve 
ne 17 COSAT! CODES 18 SUBJECT TERMS (Continue on reverse of necessan and ident, by block number! a 
Klectron density profiles Airglow 7 
lonosphere Theoretical model a 
The — 
19 ABSTRACT Continue on reverse if necessars and identity by Blok Mander = 
Since current empirical mo fels specitving tow -latitude electron densitv profiles ~ 
severely underestimate the favtinie olasnia density scale heieht and total electron content ws 
(TREO), a Semi enipirical, Low -latitude, Tomosvheric Madel (SELIM) was developed which ~. 
is notconly more realistic but is also comooutationallys fast, Electron density profiles 
(180 to 1800 km) are theoretically calo ulate d as a funtion of latitude (every 2°? between . 
24° Nand 24°S dip latitude) and local ties cesers halt hour over 24 hours LP) by solving Aa 
‘he time dependent plasmia continuity equation,  Vssuring a Chapman-like profile, sets of 
ry coeflicjents are then generate d whi ho oeprodace these individual orofiles., The coefficient: is 
themselves are easily stored, quis sve dand form the basis fora fast, portable, “4 
sep -emoirical cor souter code. dbs eesort des rites brreliy the mimput paranieters used * 
to theoretically caleubate: tor othes a vor dares use d to generate the coefficients, ‘ 
The SLE profiles are compare i awit ted Bent erepprpeal models for lquings, . i 
7 solarotrasxtiniur com ditions,  brnabbe ed: fersttoes. the coertietents, PRC and 6300 A sf, 
270 DiS" A BUTION AVARLABILITY OF ABN TRA “S 
UNO ASSL EG UNLIMITED = SAME ot RET ~: 
e 2b RESPONS- BLE oe a .* pa ; ou : mw 
¥ Va derson : ; ; s 
: DD FORM 1473, 83 APR PED ERE, moet [a haa Tae s 
5 at eee SSR CAT ON OCR Tet PAGE a 
) 
/ a: 
, os 
. 


tes 


~ 


eaiged 


na 


Unclassified 


SECURITY CLASSIFICATION OF THIS PAGE 


19, (Contd) 


airglow intensities are listed in tabular form for three seasons (Equinox, June solstice, 
and December solstice) and two solar cycle periods (solar maximum and solar minimum). 


Unclassifiec 
a a ee 
SECURITY CLASSIFICATION OF THIS PAGE 


ee oe 


w= 


pt PE Pa Ped 


be Ae Ay nee 


Preface 


We would like to thank Jeffrey Forbes, Bela Fejer, Leo McNamara, 
John Klobuchar and Jeffrey Baumgardner for numerous discussions and their 
helpful suggestions. We also wish to thank Celeste Gannon and Patricia Falcione 
for their help in manuscript typing and preparation, This research was supported 
in part by AFGI. Contract F19628-81-K-0051 with Boston University. Michael Hicks 
produced the graphs using the RS/1 package on the Boston University Department 
of Astronomy VAX-11/750. Thanks also to the staff of the AFGL Computer Center 


for their help. 


ACCeslO.i For 


NTIS CRAgy 
Ciic Tag 


AanOn sug 

Joctihcate, 1 
Sea 

By 

Di.t ib-tio / 


pit | 4%! 2c, 
OPcCial 


iii 


Goes 


. 
Oi 


(on 


TELA“ 


4 
re 
2? 
. 
4 
te 
ve, 
a 
= 
. 
“F 
is 
ny 
“J 
FZ 
‘ "> 
J 
ove 
Pad 
‘~ 
Contents ws 
ra 
1, INTRODUCTION 1 " 
1,1 Background 1 ve 
: 1,2 Model Calculations 3 
q 2. THE SEMI-EMPIRICAL APPROACH M1 ‘s, 
2.1 The SLIM Database 11 = 
2.2. Production Details 13 
2 3. COMPARISON WITH EMPIRICAL MODELS 15 G 
. 4. CONCLUSIONS 20 x4 
: REFERENCES 101 re 
re 
% 
* 
: Illustrations a 
* 
. e 
: 1. A Comparison Between Observed, Predicted, and Theoretically Aes 
. Calculated Total Electron Content (TEC) Values at Manila as es 
a Function of focal Time for January 1982 4 a 
. 2, fon-Neutral and Electron-Neutral Temperature Ratios as a 
, Function of Altitude 3 
3a, Vertical Box B Plasma Drift as a Function of Local Vime During P 
Z Solar Cycle Maximum and Minimum Periods for Mquinox Conditions 6 
‘ 3b, Vertical F< B Plasma Drift as a Function of Local Vimie During < 
Solar Cvele Maximum and Minimum Periods for June Solstice “S 
Conditic is 6 ‘of 
i] 
‘ 
« 
“ 
3 = 


“ ‘s e “a ‘ A x ad sf, a » * oh < 7 "* 
C ' 
"4 3 
é ? 
s, 
a 
i 
\ 
illustrations 7 
Vv 
¥ 
‘ 
3c. Vertical EF x B Plasma Drift as a Function of Local Time During 1 
Solar Cycle Maximum and Minimum Periods for December t 
Solstice Conditions 7 M 
4a. Contour Plot of Meridional Neutral Winds as a Function of Latitude 
and Local Time for Equinox Periods 10 
4b. Contour Plot of Meridional Neutral Winds as a Function of Latitude 3 = 
5 and Local Time for June Solstice Periods VW < 
. 4c. Contour Plot of Meridional Neutral Winds as a Function of Latitude * 
be and Local Time for December Solstice Periods 12 “ 
5. 4 Comparison Between SLIM -generated Electron Density Profiles s 
(o o) and the Empirical Models of Chiu (HM ------ @), and f 
: Bent (xxxx) for Equinox, Solar Cycle Maximum Conditions at . 
% 1200 LT 16 " 
‘ 6. A Comparison Between SLIM-generated Electron Density Profiles s 
{o o) and the Empirical Models of Chiu ( ------ BM), and 
Bent (xxxx) for Fquinox, Solar Cycle Maximum Conditions at 
2000 LT 17 
7. A Comparison Between SLIM-generated Electron Density Profiles » 
.” (o o) and the Empirical Models of Chiu (™ ------ @), and v4 
Bent (xxxx) for Equinox, Solar Cycle Maximum Conditions at uw 
0500 LT 18 
w 
8. & Comparison Between SLIM-generated Electron Density Profiles a 
(o o) and Profiles Measured by the Jicamarca Incoherent a 
Scatter Radar (e eee) on 5 March 1969 at 1400 LT and 2000 LT 19 = 
‘4 
if 
r) 
» 
» 
Tables 
1, Analytic Wind Equations ’ ae 
is 2. Electron Density Profiles and Modified Chapman Parameters | 
Re 
‘ 
ie’ 


gary tt 


sects 
os 


® 


_ A Semi-Empirical, Low-Latitude lonospheric Model 


lL. INTRODUCTION 


1.1 Background 


Current empirical models of the low-latitude ionospheric k-region (Llewellyn 
and Rents! Chiu;? and Rawer?) can severely underestimate the daytime plasma 
density scale-height and total electron content (TEC) when compared with actual 
observations, As illustrated in Figure I, the International Reference Ionosphere 
(IRI) described by Rawer™ vields daytime TEC values which are 50 percent lower 
than observed TEC values at Manila (dip 14.5° N) in January 1982 (McNamara‘), 
Some improvement occurs when the Bent topside model is incorporated with the IRI 
bottomside models but the predicted values are still lower than observed. Substantial 
improvement is achieved when theoretically calculated profiles (Anderson”) are 
used to predict TEC values, The reason is that vertical plasma transport by upward 
FE x B drift produces both topside and bottomside profiles which are much broader 
(thicker) than Chapman-like profiles. 

To calculate electron density profiles whenever they are needed, however, is 
prohibitively time-consuming on even the fastest computer, An alternate solution 


to this problem is to theoretically calculate electron density profiles as a function 


of latitude and local time and then generate coefficients which reproduce these profiles, 


(Received for publication @ October Lona) 


(Due to the Jarge number of references cited above, they will not be listed here, 
See References, page 101, ) 
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The coefficients themselves are easily stored, quickly retrieved and form the basis 
for a fast, portable, semi-empirical computer code which will produce realistic 
low -latitude F-region electron density profiles. The results of this technique are 


called the Semi-empirical Low-latitude Ionospheric Model (SLIM). 


OBSERVED TEC MANILA, JAN 1982 

R=137 DIP: 145°N 
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Figure 1, A Comparison Between Observed, Predicted, 
and Theoretically Calculated Total Electron Content 
(PRO) Values at Manila as a Function of Local Time 

for January P'82 


Using the techniques described by Anderson, electron densities as a function 
foaltitude, latitude and local time are calculated by solving the tinre-dependent 
: 
ton (QO J continuite equation numerically, “Lhe effects of production by phote - 


rompzation, boss through charge exchange with N, and oO, and transport by diffusion, 


y_v 8 


ate? 


neutral winds and vertical E Xx B drift are included. Because E XB drifts play such 
an important role in determining the electron density distribution (Hanson and 
Moftett® and Anderson”) and have a seasonal and solar cycle dependence (Fejer 

et al’), the calculations are carried out for equinox and solstice conditions for both 
solar cycle maximum and minimum periods. 

Once electron density profiles are obtained, they are normalized to the peak 
electron density, a eee and appropriate coefficients which reproduce the 
normalized profiles (180-1800 km) every half-hour local time and every 2° dip 
latitude are calculated, To regenerate density profiles simply requires specifying 
Nie gig internally (that is, re-use of war theoretical Nenax Values) or by externally 
using an empirical model such as CCIR,~ or observed foF2 values, and then applying 
the calculated coefficients applicable to the specified latitude and local time location, 
In the following sections, we describe briefly the various input parameters which 
were used in calculating theoretically the electron density profiles from the time- 
dependent continuity equation. The last section compares the current Chiu and Bent 
empirical profiles with SLIM profiles at specific latitudes and local times for 


equinox, solar cycle maximum conditions. 
1.2. Model Calculations 


The ion (0°) continuity equation is given by 
4 g 


ON; a 

ar ¢ Vien) ee 2, (1) 
where Ni(= N,) is the ion density; Pp. the ion production rate; Li» the ion loss rate; 
and Vi. the ion transport velocity. In the ionosphere, plasma is transported 
along geomagnetic field lines by diffusion and neutral winds and perpendicular to 
field lines primarily by E X B drift (Kendall and Pickering”). Solving Eq. (1) at 
low-latitudes is facilitated by transfurming the independent coordinates r, # and 9 
to a coordinate system parallel and perpendicular to B (see Anderson”), 


Equation (1) can be written 


6. Hanson, W.B., and Moffett, R.J. (1966) Ionization transport effects in the 
equatorial F-region, J. Geophvs, Res., wanes 


‘, Fejer, B.G., Farley, 1D. T., Woodman, R.IF., and Calderon, C. (1979) 
Dependence of equatorial F-region vertical drifts on season and solar cvcle, 
J, Geophys, Res., vec 


8. International Radio Consultative Committee (CCITR) (1°78) CCIR atlas of 
ionespheric characteristics, Rep, 340-3 Recommendations and Reports of the 
CCIR, Volt. 6, International Telecommunications Union, Geneva. 


§, Kendall, P.C., and Piekering, W.M. ¢1267) Ma 
P2-region altitudes, Planet, Space Sel. : 
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a arg VN, = P-L, -Ve IND SNe Vi (2) 
where Vv. is given by EX B/ B and Vin includes the effects of plasma diffusion and 
neutral wind, The right-hand side of the equation involves terms which are second 
order in the coordinate parallel to B, The left-hand side of Eq. (2) is the time rate 
of change of plasma density in a frame of reference which drifts with the EX B 
drift velocity. Equation (2) is solved numerically to give ion densities as a function 
of altitude, latitude and local time [see Mottett !° for a review of transformations 
and numerical solutions]. : 

The set of coefficients for the ion continuity equation is obtained from models 
of the neutral composition, ion and electron temperatures, and production, loss 
and diffusion rates as well as E X B drift and neutral wind models. Briefly, the 
models are the following: 

(1) The MSIS (Hedin!}) neutral atmospheric model is used to calculate 
Ny» O,, 
latitude and local time. An F10.7 cm flux of 180 units is chosen to represent solar 


and O densities and the neutral temperature, Ty as a function of altitude, 


cycle maximum period, and 70 units for solar minimum conditions. Average 
geomagnetic conditions were chosen by setting Be = 15 in the MSIS model, 

(2) Production and diffusion rates are similar to those used by Anderson. z 
For the photoionization rate at the top ofthe atmosphere, P .. a value of 
5.5 X to! sec! represents solar cycle maximum conditions, with a value of 
2.3x107! sec! for solar minimum conditions. The loss rate coefficient are 
adopted from Torr and Torr. ie 

(3) Plasma temperatures were chosen by adopting Part. and Doi ratios as 
a function of local time, altitude and season from the IRI model and then applying 
those ratios to the ue results obtained from MSIS. Figure 2 illustrates the ratio 


behavior for the sunspot maximum Fquinox case. 


10.) Moffett, R.J. (1979) The equatorial anomaly in the electron distribution of 
terrestrial F-region, Fund. Cosmis Phys., 4:3 13. 


11. Hedin, A.F. (1983) A revised thermospheric model based on mass 
spectrometer and incoherent scatter data: MSIS-83, J. Geophys. Res., 
88:10, 170 


f2. Torr, M.R., and Vorr, D.G. (1979) Chemistry of the thermosrhere and 
ionosphere, J. Atmos. Terr. Phys., 4797, 
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(4) The vertical ff) x B drift velocity observed by Woodman" is essential in re 
producing and maintaining the electron density distribution observed near the mag- a 
netic equator, We adopt six vertical drift models to represent, respectively, - 
: three seasons (equinox, summer and winter solstices? at solar evcle maximum anid 
. solar evele minimum periods, The different drift patterns are based primarily on a 
: the Jicamarca incoherent scatter radar observations reported by Fejer et al, ‘ We oe 
: ° have assumed the vertical drift velocity is independent of altitude. The east-west = 
component of F< B drift reported by Pejer et ail is neglected in the calculations p 
because its effect on electron density profiles is negligible (Aniessen "i, Figures 
= Sa, 3b and 3c display the diurnal variation in drift for the Equinox, June solstice : 
: and December solstice cases, for solar maximum and minimum, respectively, on 
. ‘4 
. 13, Woodman, R.F. (1970) Vertical drift velocities and east-west electric fields at ke 
. the magnetic equator, J, Geophys. Res., 75:6248,. . 
lt. Fejer, BOG... Farley, D.T., Gonzales, C.A., Woodman, RoF., andCalderon, ©. oa 
> (1981) F-region East-West drifts at Jicamarca, J. Geophys, Res,, 86:215,. os 
. _ — eal oe 
. 15.0 Anderson, DLN. (1981) Modeling the ambient low -latitude F-region ionosphere - am 
A Review, J. Atmos, Terr. Phys. , 43:753, . 
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Figure 3c. Vertical E XB Plasma Drift as a Function of Local Time 
During Solar Cycle Maximum and Minimum Periods for December 
Soistice Conditions 


(5) The geomagnetic field is represented by a dipole model whose axis is 
aligned with the earth's rotational axis. We incorporate EX B drift models appro - 
priate to the American sector where the dip equator lies 12° south of the geographic 
equator, and thus longitude effects are not addressed at this time. Further 
development of the model will include a tilted dipole geomagnetic field with declina- 
tion specified over the regions of interest. 

(6) A simple analytic formulation for the meridional neutral wind field is 
adopted which makes use of observations and theory. It is constructed on the basis 
of observed mean (DC), diurnal (24-hour), semidiurnal (12-hour), and terdiurnal 
(8-hour) harmonics at Arecibo (18°N), Millstone Hill (42° N) and St. Santin (45°), 
with Latitude variations of amplitude and phase obtained from theoretical (numerical) 


16, 17 


simulations (Forbes ), The geneval latitudinal pattern assumed is guided by 


16, Forbes, J.M. (1982a) Atmospheric Tides 1, Model description and results for 
for the solar diurnal component, J. Geophys, Res., B7;4222, 


17.) Forbes, J.™M. (1982b) Atmospheric Tides 2.0 The solar and lunar semi-diurnal 
components, J. Geophys. Res., 872524 


vr 


wiv 2 we 


Pair 


oa 
ke ae at et a! at 


s 


ee ee 


rere arate es | elt hos 


eae 


oe tee ® 


Were, 


18, Roble, R.G., Dickinson, R.E., and Ridley, E.C., (1977) Seasonal and 


the tidal calculations of Forbes! ** te 


Roble et al. es The analytic expressions, valid only for latitudes lel = 45°, are 


and the zonal mean velocities presented by 


given in Table 1 for Equinox, Dec*mber solstice and June solstice periods. 
For solstice conditions, the diurnal component, Vues shifted 10° in aa 
), The Roble et al 
mean meridional winds, Vo during solstice are shifted downwards by about 


in the direction of the solar bulge (consistent with Forbes 


40 m/sec so that the transequatorial amplitudes are of the order 50 in/sec and 
the winter hemisphere reversal occurs near Millstone Hill. 

The solsticial modifications to the mean and diurnal components adequately .- 
account for the main features of the observed seasonal-latitudinal dependence of the 
meridional wind for our purposes, Seasonal variations in the semi-diurnal and 
terdiurnal components are neglected, Figures 4a, 4b, and 4c give contour plots of 
the overall wind model spanning the simulation region where positive values 
represent north-to-south winds, Finally, it should be noted that while EUV varies 
markedly over a solar cycle, the associated variations in ion drag suppress solar 
cycle variations in the tidal winds, although neutral temperatures vary consider- 
ably, Therefore, we have assumed the current wind model parameterization to 
be independent of solar activity. 


and solar cycle variations of the zonal mean circulation in the 
thermosphere, J, Geophys. Res., Besotee 
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Table 1. Analytic Wind Equations 


6 is latitude in degrees 
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Table 1. Analytic Wind Equations (Contd) 
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Figure 4a, Contour Plot of Meridional Neutral Wind as a kunction 


of Latitude and Local Time for Equinox Periods, Positive values 


represent north-to-south winds 
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Figure 4b, Contour Plot of Meridional Neutral Wind as a Function 
of Latitude and Local Time for June Solstice Periods, Positive 
values represent north-to-south winds 


2. THE SEMI-EMPIRICAL APPROACH 


2.1 The SLIM Database 


As described in Section 1. 1, the theoretical model is used to produce six pre- 
liminary databases. Three of these runs are for sunspot minimum (FF 10.7 = 70, 
Ap 15) during the June solstice, December solstice, and Equinox seasons. The 
remaining three are for sunspot maximum (F 10.7 180, Ap - 15) during the same 
seasons. Each record of data consists of profile parameters stored for a fixed 
geographic: (geomagnetic) point and fixed time, Fach database covers 25 latitudes 
(24° 5 - 24° N, every 2?) and 48 times (0000h - 2330h, every 30 min). A subset of 
the SLIM data set, using a more coarse resolution (NCH) every 4° in latitude and 
every hour of local tims) is presented in tabular format for easy access anid 


reference (Table 2), 
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Figure 4c. Contour Plot of Meridional Neutral Wind as a Function 
of Latitude and Local Time for December Solstice Periods. 
Positive values represent north-to-south winds 


The core of the SLIM approach is to identify a relatively simple N,(h) function 
that will have a small number of coefficients and yet provide a good representation 
of the profiles calculated from the theoretical model. As guiding principles, it was 
decided that; (1) maximum accuracy was required near the I-region peal: (in both 
density and height), and (2) while some errors in No could be accepted well below 
or above hee the overall profile (as represented by the total electron content 
(ThC) parameter) had to reproduce the theoretical results to within 5 percent or less. 
The Chapman profiles developed in the earliest days of ionospheric physics 
offer a simple way to explain the vertical structure of plasma in the upper atmosphere, 
They are based on relatively simple physics and thus cannot be used to represent 
the entire Nth) profiles that emerge from the more sophisticated simulation tech - 
niques in use today, As simple mathematical formulas, however, they are well- 


suited for representing portions of the ionospheric profile. Specifically, while a 
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single Chapman profile cannot fit entire N,(h) profile, separate Chapman formulas 
can be used to reproduce the topside and bottomside components of the overall 
profiles with good results, it was decided, therefore, to define a ''generalized" 

or "modified" Chapman formula with coefficients that vary for each topside and 
bottomside component of the N,(h) profile. The function used to fit the theoretically - 
terived profiles was thus chosen to be; 

Zz 


N,(h) = Ne vax exp (c(1 - z-e 
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(h - bax’ 
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In using Eq. (3), six parameters (or coefficients) are needed to reproduce the 
N,(h) profile: re Hes soe ue Aup Co: Ajo The advantages of this method 
are: (1) the fit achieves 100 percent accuracy at Nag and Ms a since they are 
themselves input coefficients; (2) the topside and bottomside Chapman-type 
functions guarantee continuity of both density and its first derivative across the 
fit boundary; and (3) with a fit technique weighted by density values (see below), 
the c's and A's chosen emphasize maximum accuracy in regions near the peak. 

Included on each record of the database are two other numbers available for 
use and/or inspection. These are; (1) the total electron content (TEC) obtained 
by integrating the SLIM profile from 180 to 1800 km; and (2) 6300 A emission 
calculated from os recombination chemistry only (and therefore, appropriate for 
post-sunset airglow), 


2.2 Production Details 


There are four steps and one version of the database is produced after each 
step. Initially, the theoretical model LOLAT is used to produce Database 1 (DB1). 
DBI simply contains the electron density and accompanying heights. The model 
solves the o* continuity equation, fieldline~by-fieldline, starting at the highest 
altitude, The order of the density and height arrays are latitude, time, and 
fieldline. The computer file is sequential. 

Database 1 is not organized into actual electron density profiles. It is the job 
of program SEQDIR to sort the densities by height, time, and latitude, In addition, 
SEQDIR creates Database 2 (DB2) as a direct access file. A direct access file is 
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one in which each record of the file can be individually accessed without reading 
any other record, In DB2, each record contains a profile of electron densities 
and a height array for a particular time and latitude. 

Iatabase 2 has the disadvantage that the electron densities are given at 
unevenly spaced height intervals. EVEN is the program used to interpolate between 
the given values. The resulting profile records contain values at every 20 km from 
180 to 1000 km and every 100 km from 1000 to 1800, The interpolation is done with 
the National Center for Atmospheric Research (NCAR) family of cubic spline inter- 
polation, derivative and integration routines called CURV. The derivative routine 
is used to determine the peak electron density by identifying heights where the 
derivative changes sign. Arrays of the densities at these heights are scanned for 


the maximum value, which is identified as Nun The density profile is stored 


at predetermined heights so that the height poet is no longer necded, making 
LLatabase 3 (IDB3) much smaller than DB2, TEC from 180 to 1800 km is also 
calculated using the CURV routines and included in each DB3 record, 

Database 4 is the result of program CHAP fitting the profiles in DB3 to the 
generalized Chapman formula [Eq. (3)}]. Nak and Ris ge are known from DB3. 
Combinations of values of c and A are tested until the sum of the square of the 
differences, 
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is minimized, The range of values for c and A in the first loop is 


4 = +2.0, Ace a: A * 10,0 + 200.0 km, AA = 10.0km. 


If the minimum occurs at the endpoint of either range, the range is extended 
and the cycle is repeated with the same Ac and SA, Once an initial c and A are 


selected (c', A'), a new cvcle is started where 


1 ; 1 Eien? val 
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As (A' - 10.0) > (A' + 10.0), OA= 1.0. 


To save space, the C quoted in the SLIM tables is actually 


C= 1000 * c. 
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The final c and A values are added to DB4 and TEC is recalculated, he 
The end product is the SLIM database (DBSLIM), the ensemble of the DB4's for 


the six conditions of solar cycle and season, 


3. COMPARISON WITH EMPIRICAL MODELS 
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Figure 5 compares the Semi-empirical, Low-latitude, Ilonospheric Model (SLIM) 
with the Chiu and Bent models for Equinox, solar maximum conditions at 1400 LT. 
: The daytime upward EX B drift pictured in Figure 3a causes the broad F-region pro- 
file. Respective topside scale heights at 700 km are 540 km, 150 km and 115 km. 
Integrating the electron density profiles with altitude give total electron contents (TEC) fn 
of 97, 48, and 61 (x 101° ake), respectively. The bottom portion of Figure 5 


shows the comparison at 16°S dip latitude at the crest of the equatorial anomaly. 


= a ys “ * 


« 


Although the SLIM topside scale height is not as large as at the magnetic equator, it 


is still larger than those given by Chiu or Bent. In addition, the Nan value of 
3,5 X 10° Sone is significantly larger than the value at the magnetic equator due to 


the combination of upward and poleward E X B drift and downward diffusion, Note 
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that neither the Chiu or Bent profiles show any significant change with latitude, 


The most dramatic differences between SLIM and the empirical models occurs 


just after sunset around 2000 LT and pictured in Figure 6. At the magnetic equator 


{top portion), the post-sunset enhancement in upward drift lifts the F-layer to 600 km 
‘ompared to eee values ot 450 km and 350 km for Bent and Chiu models, HEspechivelys 
Such differences in Pe oa would cause significant differences in the estimated 6300 A 
airglow intensity and the ion-neutral drag coefficient in the neutral momentum equation. o 
At 16°S dip latitude (bottom portion of Figure 6) the same upward drift causes an 
- enhancement in N . The calculated value is 5 X 108 elfem® compared with a 
* 2x 108 eiverne for Chiu and 8 X 10” ethom for the Bent model. Respective TEC 


p values (X 10}? 


el/em?) are 139, 48 and 19. Finally, Figure 7 gives a comparison 
. between SLIM, Chiu and Bent at 0500 LT. The calculated profiles at the magnetic 
equator and 16’S dip latitude are much more Chapman-like and in reasonable ayree- a 


ment with the empirical models in both N h and topside scale height. The cal- 


re j max’ max 12 3 
culated value of ‘TEC at 16°S dip latitude (bottom portion) is 5 X 10°" el/em* which 
: compares favorably with the Bent value of 3.5 X 10! and a Chiu value of 9X 10!” elven. . 
Figure 8 compares electron density profiles measured on 5 March 1969, using , 

the Jicamarca incoherent scatter radar facility with calculated density profiles at = 

. 1400 LV and 2000 LT. Yhe comparison demonstrates that both the large davtime ‘ 
; plasma scale height and the calculated value of Lays at 2000 LT are in reasonable ‘ 
agreement with observed values, Again, it should be emphasized thai the profiles : 


are sensitive to the vertical Ex B drift velocity whose diurnal variation mav have 


significant day-to-day changes, 
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Figure 8 A Comparison Between SiLU]-generated Mlectron 
Density Profiles (0 ——o) and Profiles Measured by the 
Jicamarca Incoherent Scatter Radar (ee ee) ond March 1969 
it 1400 LT (top portion) and 2000 LV (bottom portion) 


4. CONCLUSIONS 


A very brief description of the generation of a Semi-empirical, Low-latitude, 
Jonospheric Model has been presented in which theoretically calculated electron 


density profiles are approximated by a Modified Chapman function given by Iq. (4). 


t fs 2 = -Z 
N(h) 7 Ne exp ([c(l-z-e -)] (4) 
where 
ba h - Amax 
—- 


and N and h represent the peak electron density and the altitude of the F2 
m max heh 


peak, meee e are Above the peak altitude, aye = = real and c= Sup and below 

2 max ee Bes oe se : , ; ; 
peak, 257 mo rae and ¢ Clot The six coefficients which describe each profile 
are IN sed cos eee Aup pe Ato and C16 Appropriate coefficients which regenerate 


electron density profiles (180-1800 km), every half-hour local time and every 

2° dip latitude between 24° N and 24° dip latitude have been calculated. Six cases 
have been generated, both solstice periods and one equinox period for solar cycle 
minimum and solar maximum conditions. 

A subroutine which is able to quickly generate realistic electron density profiles 
has a number of applications. By providing diurnal N, (h) and TEC values at any 
specified dip latitude or slant path, SLIM could be used to examine propagation 
effects determined by ionospheric structure. SLIM has several possible applica- 
tions related to so-called “active experiments" where ambient ionospheric condi- 
tions need to be specified for pre-experiment planning and post-event analysis. In 
addition, such a model would yield flux-tube integrated Pedersen conductivity and 
electron content values which are important in determining and predicting low - 
latitude, instability growth rates. Finally, the program would be capable of 
supplying airglow intensities both in the vertical and at any slant viewing angle 


desired, 
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Table 2, 


Electron Density Profiles and Chapman Parameters 


Each page of Table 2 provides the following values for a 


specified latitude, season and solar cycle: 


1, Electron densities Mo fen) from 180 to 1800 km. 


‘ Y 3 es 3 
2. Mare (10° /cem*~) and Dias: (km). 


3. The four Chapman parameters Bs (km), Sa 
A, (km) and C 
lo 


lo’ 


2 D) 
4, Total electron content - TEC (10!2/ em?) and 


630 nm o* recombination airglow (Rayleighs), 


The above values are listed every hour local time and 


cover the latitude range from 24° N dip latitude (+24°) to 


24°S dip latitude (-24°) in 4°-increments. The table is 


grouped by solar cycle (maximum and minimum) and within 


each solar cycle by season (Equinox, June Solstice and 
December Solstice), Remember that C), = 1000 * on 


and C, = 1000% ¢ 
lo 


, where ap and S69 are defined by Eq. (3), 
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